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The possible heavy fermion superconductivity is carefully reexamined in the two-dimensional Kondo lattice 
model with an antiferromagnetic Heisenberg superexchange between local magnetic moments. In order to 
establish an effective mean field theory in the limit of the paramagnetic heavy Fermi liquid and near the half- 
filling case, we find that the spinon singlet pairing from the local antiferromagnetic short-range correlations can 
reduce the ground state energy substantially. In the presence of the Kondo screening effect, the Cooper pairs 
between the conduction electrons is induced. Depending on the ratio of the Heisenberg and the Kondo exchange 
couplings, the resulting superconducting state is characterized by either a d-wave nodal or d-wave nodeless state, 
and a continuous phase transition exists between these two states. These results are related to some quasi-two 
dimensional heary fermion superconductors. 
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I. INTRODUCTION 

Heavy fermion materials have been playing a particular im- 
portant role in our understanding of strongly correlated elec- 
tron systems[i, 2], and the Kondo lattice model is believed to 
capture the basic physics of heavy fermion systems[3]. The 
model describes a lattice of local spin- 1/2 magnetic moments 
coupled antiferromagnetically to a single band of conduction 
electrons. When the number of conduction electrons is less 
than the number of the local magnetic moments, the coher- 
ent superposition of individual Kondo screening clouds gives 
rise to the huge mass enhancement of the quasiparticles, and 
the resulting metallic state is characterized by a large Fermi 
surface with the Luttinger volume containing both conduction 
electrons and local moments [ , ] . Competing with the Kondo 
singlet formation, the local magnetic moments indirectly in- 
teract with each other via magnetic polarization of the conduc- 
tion electrons - the Ruderman-Kittel-Kasuya- Yosida (RKKY) 
interaction. Such an interaction dominates at low values of the 
Kondo exchange coupling and is the driving force for the an- 
tiferromagnetic (AFM) long-range order near the half-filling 
of the conduction electrons [6-8]. 

In addition, there have been growing evidence that the an- 
tiferromagnetism is also intimate with superconductivity in 
some typical heavy fermion compounds, such as CeCu2Si2 
(Ref.[9]) and CeRhlns (Ref.[ ]). So far various mechanisms 
for heavy fermion superconductivity have been studied, in- 
cluding paramagnon exchange[l 1, 12], conventional phonon- 
mediated[ ], and Kondo-boson-mediated pairings [ - ]. 
So far the pairing mechanism of the heavy quasiparticles is 
still under investigation theoretically. However, various exper- 
imental evidences in some heavy fermion materials strongly 
support the d-wave pairing superconductivity [ ] . 

It is well established that the large- fermionic approach 
can be used to treat the Kondo lattice model very efficiently, 
leading to a stable paramagnetic heavy Fermi liquid state [ , 
15, 18]. Qualitatively, the RKKY interaction promotes sin- 
glet formation among the local magnetic moments, reducing 



the tendency of singlet formation between the local moments 
and conduction electrons. To simplify the RKKY interaction, 
one can explicitly introduce the local AFM Heisenberg su- 
perexchange Jh among the local moments to the Kondo lat- 
tice system[ 19-25]. Thus, the paramagnetic heavy Fermi liq- 
uid state {Jk > Jh) provides a good starting point for theo- 
retical considerations, because a number of instabilities can 
be further discussed, including the AFM ordered state and 
unconventional heavy fermion superconductivity. Recently, 
some numerical evidences on robust d-wave pairing correla- 
tions have been given in small size Kondo-Heisenberg lattice 
system[^^^]. 

In this paper, we try to establish an effective mean field 
(MF) theory to the Kondo-Heisenberg lattice model on a two- 
dimensional square lattice in the limit of Jk ^ Jh- After 
carefully examining three different MF treatments for the lo- 
cal AFM Heisenberg superexchange, we find that the spinon 
or f-fermion pairings can substantially save the ground state 
energy in the presence of the Kondo screening. Near the 
half-filling, the Cooper pairs between the conduction elec- 
trons can be induced via the Kondo screening effects. We 
further find that the resulting superconducting pairing func- 
tion has a d-wave symmetry. Whether there are nodes or not 
depends on the ratio of the local AFM Heisenberg and the 
Kondo exchange couplings x = Jh/ Jk- There is a contin- 
uous phase transition between the nodal and nodeless super- 
conducting phases. Before the magnetic phase transition, the 
Kondo screening order parameter is strongly suppressed, and 
the present mean field theory is no longer valid. 



II. MEAN FIELD TREATMENT 

The model Hamiltonian of the Kondo-Heisenberg lattice 
model is given by: 

H = Y^ 6k4^Cka ^ JR^Si'^i^ Jn^Si' S^, (1) 
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where cj.^ creates a conduction electron on an extended or- 
bital with wave vector k and z-component of spin a =t, 
The spin- 1/2 operators of the local magnetic moments have 
the fermionic representation = | Xlcrcr' flayer a' fia' where 
r is the Pauli matrices. There are two local constraints: 
ficjfi^y = 1 fi^fii = 0- The former local constraint 
restricts any charge fluctuations, the f-fermions just describe 
the spin degrees of freedom of the local moments, and we will 
refer them to spinons. The latter constraint is imposed by the 
spin SU(2) symmetry and the extended s-wave spinon pairing 
between spinons can be excluded. So the paramagnetic Fermi 
liquid limit (Jk ^ Jh) provides us a good starting point of 
theoretical considerations. 

It is first noticed that the Kondo spin exchange can be ex- 
pressed as the singlet pairing attraction between the spinon 
holes and conduction electrons up to a chemical potential shift 



(2) 



Then the Kondo screening order parameters can be introduced 
as 



(3) 



For the non-magnetic states, the local AFM Heisenberg su- 
perexchange can be expressed in terms of either the spinon 
hopping or singlet pairing operators 
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However, most of previous theoretical studies on the Kondo- 
Heisenberg lattice model [ - ], only introduced the spinon 
hopping order parameter. Actually, Coleman and Andrei [19] 
had emphasized that the local SU(2) gauge invariance of the 
local Heisenberg spin operator generally requires the consid- 
eration of both MF order parameters. Moreover, recent ad- 
vances in this area have been made by using symplectic repre- 
sentation of the local magnetic spins [ ] . Here we first intro- 
duce two MF valence bond order parameters to characterize 
the short-range AFM correlations between the local moments 

Xii = -(/^/it + /i/.i), A., = -if^l - /^/|^). (5) 

To avoid the incidental degeneracy of the conduction 
electron band on a square lattice, we choose ek = 
—2t (cos kx + cos ky) + At' cos kx cos ky — ji, where t and 
t' are the first and second nearest neighbor hoping matrix 
elements, while ji is the chemical potential, which should 
be determined self-consistently by the density of conduc- 
tion electrons ric. When the local AFM correlation hop- 
ping parameter x is simply chosen as a uniform parame- 
ter, the spinons form a very narrow band with a dispersion 
Xk = -^i^x(cos kx + cos ky) + A, where A is the Lagrangian 
multiplier to impose the local constraint flafio- = 1 on 
average. For the short-range AFM spinon singlet pairing or- 
der parameter A^j, the local pairing constraint = 



excludes the extended s-wave pairing, and then d-wave sym- 
metric pairings lead to a lower ground state energy [2 8], cor- 
responding to Ai,i+e, = -/^i^i+ey = Aq. 

Then the MF model Hamiltonian in momentum space can 
be written as 
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where Ak = Ao(cosA:a; — cos ky) and eo = 
Jh^o + JhX^ ~ When a Nambu spinor is defined by 

V^k = ( 4t ^-H fit f-^^ )' express the MF 

model Hamiltonian in a matrix form 
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Diagonalizing this MF model Hamiltonian, we obtain two 
quasiparticle energy bands 



k2' 



(8) 



where 



^k2 = \ (£kXk - + (Sk JffAk)2. (9) 



Due to the particle-hole symmetry of the superconducting 
quasiparticles, all the negative energy states are filled up at 
zero temperature, and the ground state energy per site is thus 
given by 



JkV' 



J^(Ag + x'), (10) 



where E^^ = E:^ ^E^ = y/2E^[T2E^. The saddle point 
equations for the MF order parameters V, x, Aq and A can be 
determined by minimizing the ground state energy: 



dEa dEa dEa dEa 



dV 



dx 



dX dAo 



0. 



(11) 



The chemical potential /ic is determined by the relation Uc = 
— . Therefore, the self-consistent equations at zero tem- 
perature are given by. 
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In order to simplify the present treatments, we can have 
three different MF schemes: (1) x = and Aq 7^ 0; (2) x 7^ 
and Ao 7^ 0; (3) X 7^ and Aq = 0. Without losing general- 
ity, we choose t' jt = 0.3, Uc = 0.8 and fixing Jj^/t = 2.0, 
we have carefully solved the self-consistent equations and 
compared the ground state energy densities for the three dif- 
ferent types of MF schemes. The numerical results are not 
sensitive to the parameters chosen. The numerical results are 
displayed in Fig. 1. It is clear that, for < J/^/t < 1.3, 
the first type of MF ansatz has a comparatively lower ground 
state energy. As the local AFM coupling J/^/t increases, the 
ground state energy is almost unchanges, and then decreases 
gradually as Jh/'^ is further increased. However, the ground 
state energies for two other MF schemes grow up in the begin- 
ning and then blend down after reaching their maximal points. 
Only when Jh/^ > 1.34, the ground state energy of the sec- 
ond MF ansatz starts to become lower than that of the first 
MF ansatz. In the present paper, we will confine our follow- 
ing discussions to the small local AFM Heisenberg exchange 
coupling Jn/t < 1.34, so only the first type of MF scheme 
(X = and Aq 7^ 0) will be used. 
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FIG. 1: (Color online) The ground state energies of a two- 
dimensional Kondo-Heisenberg lattice model under three different 
types of MF approximation schemes for a given value of Jx/t = 
2.0. The other parameters are t' /t = 0.3 and Uc = 0.8. 



III. D-WAVE SUPERCONDUCTING PHASE 

When there is only one short-range AFM correlation order 
parameters Aq 7^ 0, the MF model Hamiltonian is simplified 
by 

kcr kcr 
kcr k 

+ N{^ + JHAl-Xy (13) 

where the spinons have a flat band. Due to the hybridiza- 
tion effect, the heavy quasiparticle bands form. On the other 
hand, the spinon singlet pairings imply a special non-magnetic 



state with short-range AFM correlations among the local mag- 
netic moments. Although there are no any direct attractions 
among the conduction electrons, the spinon singlet pairings 
provide an indirect glue for the formation of the Cooper pairs 
via the Kondo screening/hybridizing effect. So the resulting 
MF ground state represents a heavy fermion superconducting 
state. 



On a two-dimensional square lattice, the spinon singlet 
pairing function has a d-wave symmetry [ ]. Whether there 
are nodes or not mainly depends on the Lagrangian multiplier 
A and hybridization strength V. From the simplified quasipar- 
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FIG. 2: (Color online) The lower branch of the superconducting 
quasiparticle band struture for different local AFM Heisenberg cou- 
pling strength with t' jt = 0.3, ric = 0.8 and Jx/t = 2.0. 



tide bands 
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(£kA-^^)2 + (ekJifAk)2, 



(14) 



we note that nodes appear in the lower band when £^k,2 = 

0, which requires the condition £k = diagonal 
direction. 

Then the numerical mean field calculations are carefully 
performed with the parameters t^t = 0.3, ric = 0.8 and 
Jx/t = 2.0. The lower branch of the quasiparticle band 
in the first Brillouin zone along the momentum direc- 
tion (0, 0)^(7r/2, 7r/2)^(7r, 7r)^(7r, 7r/2)^(7r, 0) are plotted 
in Fig. 2 for different local AFM coupling strength Jn/t = 
0.5, 0.8, 1.192, and 1.3, respectively. We first notice that the 
quasiparticle states near the Fermi energy have a very small 
dispersion, reflecting the quasiparticle mass enhancement in 
the superconducting state. In Fig. 2a and Fig. 2b, a node can 
be clearly seen between the momentum (7r/2, 7r/2) and (tt, tt) 
in the diagonal direction. Secondly, the position of the node 
changes as increasing the local AFM coupling strength. When 
the nodal position is shifted to (tt, tt), further increasing Jn/t 
makes a small energy gap open. The critical value can be de- 
termined as JH/t = 1.192. When 1.192 < Jn/t < 1.34, 
the resulting superconducting state has a full superconducting 
gap. 

Moreover, the Kondo screening parameter V and the AFM 
spinon pairing parameter Aq have also been calculated, and 
the numerical results are shown in Fig. 3. It can be seen that, as 
Jn/t starts to increase, the Kondo screening parameter is al- 
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FIG. 3: (Color online) The effective mean field order parameters of 
the Kondo screening V, the f-fermion/spinon pairing Ao, and the 
Cooper pairing Asc for a given value of Jk /t = 2.0. 



most unchanged in the beginning and then becomes gradually 
decreased when Jn/t > 0.80. On the contrary, the spinon 
pairing order parameter Aq is extremely small in the small 
limit of Jn/t. Only after Jh /t > 0.20, Aq starts to grow up 
quickly. Actually, we have also plotted a superconducting or- 
der parameter Asc defined by the conducting electrons, whose 
definition will be given later. It should be pointed out that the 
quantum phase transition mentioned above does not manifest 
in the MF order parameters. 

As the spinons are paired up in the real lattice space, it 
is more interesting to display the spinon pairing distribution 
in the momentum space. Using the method of the equation 
motion, we can easily derive the double-time retarded Green 
function for the spinon pairs. 



kil/kt//^+^o+ 



Ji^Ak(c^^ -4) 



(15) 

from which the expectation value of the spinon pairing func- 
tion can be evaluated according to the spectral theorem of the 
Green's function. At zero temperature, we obtain 



f\ ft 

^kt-^-ki 



JhA^ 



2^2 (^ki+^k2) L 



1 



E\^2 



(16) 



which is displayed in Fig. 4. When the local AFM coupling 
strength Jj^/t = 0.5, we observe in Fig.4a that the spinon 
pairing function reflects the d-wave symmetry around the 
Fermi surface of heavy quasiparticles, which is a hole-like 
circle around the corner of the first Brillouin zone. A node 
is clearly seen between the momentum (7r/2, 7r/2) and (tt, tt) 
in the diagonal direction. The positive pairing magnitude and 
the negative pairing magnitude are separated by the nodes. As 
increasing of Ji^^/t, the spinon pairing region is outstretched 
and the position of the node is shifted toward to (tt, tt), as 
shown in Fig.4b. When Jj^/t = 1.192, the node starts to 
disappear, and a small energy gap opens at (tt, tt). 
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FIG. 4: (Color online) The spinon singlet pairing distribution in the 
first Brillouin zone as increasing local AFM coupling strength for a 
given value of Jx/t = 2.O.. 



However, the Cooper pairing function of the conduction 
electrons is more important for the heavy fermion supercon- 
ducting state. Similarly, we can easily derive the double-time 
retarded Green function for the conduction electron pairs 



(17) 

and the superconducting pairing function can be defined from 
the expectation value of the conduction electron pairs. At zero 
temperature, we have 



8£^k2\/2 (£^ki + E\^2) 



(18) 



which has been displayed in Fig. 5. When Jh /t = 0.5,we can 
clearly observe that the positive pairing is distributed around 
the upper part of the Fermi surface, while the negative pair- 
ing is around the lower part of the Fermi surface. Nodes 
sit in the diagonal directions of the Fermi surface. In con- 
trast to the spinon pairing distribution, there are additional 
features: a positive electron pairing magnitude is centered 
around the point (0, tt) and a negative electron pairing mag- 
nitude is around the point (tt, 0). As the local AFM coupling 
Jh/^ > 1.192, an energy gap opens near the nodes of the 
Fermi surface, and there are only small pairing distribution on 
other parts of the Fermi surface. However, the Cooper pairing 
distributed around momentum (0, tt) and (tt, 0) have strong 
magnitudes. The resulting superconductivity still exhibits the 
d-wave symmetry, i.e., the d-wave nodeless superconducting 
state. 

Moreover, a real superconducting order parameter can be 
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FIG. 5: (Color online) The superconducting pairing distribution 
function in the first Brillouin zone for the local AFM coupling 
strength Jn/t = 0.5 and 1.192 with a given value of Jx/t = 2.0. 
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FIG. 6: (Color online) The ground state energy density (a) and its 
first order derivative (b) as a function of Jh /t for a given value of 
Jx/t = 2. The arrow indicates the critical point of quantum phase 
transition. 



defined by 

^sc = (cosA:^ - cos ky) {cl^c^_^^ + C-k^Ckt), 

k 

(19) 

which has been displayed in Fig. 3. Since it is generated by the 
spinon singlet pairings, the superconducting order parameter 
has a smaller magnitude compared to the spinon pairing order 
parameter Aq. In order to reveal the quantum phase transition 
between two superconducting phases, we have to calculate the 
first order derivative of the ground state energy. For a given 
value of Jx/t = 2.0, the numerical results are exhibited in 
Fig. 6. At Jn/t = 1.192, the ground state energy and its first 
derivative are continuous but the second derivative is not. This 
corresponds the critical point of a second order quantum phase 
transition. 



IV. DISCUSSION AND CONCLUSION 

So far we have presented the ground state properties of the 
effective MF theory for the two-dimensional Kondo lattice 
model with the local AFM Heisenberg exchange coupling be- 
tween the localized magnetic moments. We would like to em- 
phasize that, it is the local AFM short-range interaction that 
induces the AFM and superconducting long-range ordering 
states. When Jk ^ Jh, the local AFM Heisenberg exchange 
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coupling among the local moments is negligible and the sta- 
ble paramagnetic heavy Fermi liquid state is resulted from the 
coherent superposition of individual Kondo screening clouds. 
When Jh is increased but still less than J^, the local AFM 
Heisenberg exchange coupling has to be taken into account, 
and the spinon singlet pairings can reduce the ground state 
energy significantly. The Cooper singlet pairs among the con- 
duction electrons can be induced via the Kondo screening ef- 
fect, leading to the heavy fermion superconductivity. This also 
belongs to the category of the unconventional superconductiv- 
ity mediated by the short-range AFM correlations. However, 
in the present MF theory, we do not find a finite critical value 
of JH/t to separate the heavy fermion liquid and the super- 
conducting phases. This may be caused by the approximation 
used in the effective MF theory. 

The resulting superconducting state exhibits the d-wave 
symmetry. Whether the Cooper pairing function has nodes 
or not depends on the ratio of the AFM Heisenberg and 
Kondo exchange couplings. As the coupling ratio increases, 
the nodal position is shifted outward along the direction of 
(7r/2, 7r/2)^(7r, tt). When the nodal position reaches the end 



point of the Brillouin zone, a quantum phase transition occurs, 
and a full superconducting gap is opened at the Fermi energy. 
As further increasing of Jh/ Jk, the present MF theory is no 
longer valid, because the long-range AFM correlations have 
to be taken into account. 

In conclusion, a possible mechanism of heavy fermion su- 
perconductivity with d-wave symmetry is carefully investi- 
gated in the two-dimensional Kondo-Heisenberg lattice model 
from the limit of the paramagnetic heavy Fermi liquid Jk ^ 
Jh- The resulting d-wave superconducting properties can be 
related to some heavy fermion superconductors with the sim- 
ilar structure like CeCoIns, where thermal conductivity mea- 
surements strongly support a superconducting gap with nodes 
along the diagonal directions in the Brillouin zone[29]. Fur- 
ther theoretical work including the estimation of the fluctua- 
tions around the MF solution or the coexistent phase of both 
the AFM order and superconductivity would be necessary to 
be considered. 
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